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We present one-loop results for the forward twist-two matrix elements relevant to the unpolarised, 
helicity and transversity baryon structure functions, in partially-quenched (Nf = 2 and Nf = 2 + 1) 
heavy baryon chiral perturbation theory. The full-QCD limit can be straightforwardly obtained 
from these results and we also consider SU(2|2) quenched QCD. Our calculations are performed in 
finite volume as well as in infinite volume. We discuss features of lattice simulations and investigate 
finite volume effects in detail. We find that volume effects are not negligible, typically around 5-10% 
in current partially-quenched and full QCD calculations, and are possibly larger in quenched QCD. 
Extensions to the off-forward matrix elements and potential difficulties that occur there are also 
discussed. 



I. INTRODUCTION 



The quark and gluon sub-structure of hadrons has been probed for many years in high energy scattering experiments. 
Much of the information that has been gleaned is encoded in the parton distribution functions (PDFs) that describe 
the longitudinal momentum distributions of quarks and gluons within hadrons. Cross-sections for deep inelastic 
scattering, for example, have been shown to factorise into short and long distance contributions 0,0- The short 
■ distance pieces (Wilson coefficients) are perturbatively calculable whilst the long range effects are expressed in terms 
of the PDFs. The utility of such PDFs is that they are universal; the same set of PDFs appear in deep-inelastic 
t~ — \ scattering, Drell-Yan processes and heavy vector boson production. Whilst PDFs are scale dependent, once they are 
known at one scale they can be calculated at higher scales via the DGLAP [3j evolution equations. A number of 
groups 0, || have exploited the universality of PDFs and their known scale-dependence by performing global analyses 
of experimental data, thereby providing convenient parameterisations of the PDFs. Such parameterisations have 
proven very useful in testing perturbative QCD in high energy processes and constraining new physics, but nothing 
, is learnt about the non-perturbative origins of the PDFs. 

Whilst experiments continually increase our knowledge of PDFs, there is much that is still unknown. Recent results 
IE 111 have shown that u(x) ^ d(x), however other simple qualitative questions such as whether Au(x) = Ad(x) 
• or s(x) = s(x) remain unanswered. Even for the unpolarised valence quark distributions, information is scarce at 
i' large x and there is no experimental information about the transversity distributions. Consequently, any insight that 
can be gained directly from QCD would be very useful. Since the PDFs encode the soft, hadronic scale physics of 
QCD bound states, perturbative QCD is of little use. One can turn to models to suggest the qualitatively important 
features of PDFs (for example u(x) ^ d(x) was predicted on the basis of the pion cloud |9|), but to make concrete 
predictions with systematically improvable errors one must solve QCD non-perturbatively. Currently this means one 
must use lattice QCD. 

In lattice QCD, one discretises space-time and uses Monte-Carlo techniques to evaluate the functional integrals over 
the quark and gluons fields, necessarily making a Wick rotation to Euclidean space in the process. However, deep- 
inelastic scattering and related processes are dominated by distances that are light-like, and as such are inaccessible 
in Euclidean space calculations. The way around this difficulty is provided by the operator product expansion (OPE) 
which relates matrix elements of certain local operators to Mellin moments of the various quark and gluon distributions 
(defined below). For quark distributions, the twist-two (twist = dimension — spin) operators that arise from the OPE 
of the bilocal light-cone operators in Nf — 2 QCD are 



QCB n (A) 
QCDfi(A) 

^■■■fn 

QCD qT;{A) 



V'Tf/xo^Mi ■ ■ • D n n } T A?P - traces 



(1) 



^7(^0 75-0^1 • ■ • D^taiP - traces , (2) 



tyvnaiiiilsD^ . . . D^ ln+l} T A ip - traces 



(3) 



where ta is an isospin matrix (to = 1, 1X2,3 are the Pauli matrices), {. . .} indicates symmetrisation of indices, the 

gauge covariant derivative = \(D^ — D^), and 'traces' are subtracted in order that the operator transforms 
irreducibly under the Lorentz group. Additional twist-two operators can be built exclusively from gluon fields and 
towers of higher twist operators can also be constructed, but we shall not consider these here. 
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The above operators are related to the spin averaged, longitudinally polarised, and transversely polarised quark 
distributions . We first define the Mellin moments of the quark distributions q = q' +q, Aq = q' — q^- and Sq = q T —q 1 - 
(where qTU) corresponds to quarks with helicity aligned (anti-aligned) with that of a longitudinally polarized target, 
and q T ( ± ' > corresponds to quarks with spin aligned (anti-aligned) with that of a transversely polarized target) for 
flavour q as 

(x n ) q = [ dxx n [q(x)-(-l) n q(x)] , 
Jo 

(x n ) Aq = / dx x n [Aq(x) + (-l) n Aq(x)} , (4) 
Jo 

(x n ) Sq = [ dx X n [Sq(x) - (-l) n Sq(x)} . 
JO 

These moments are then related to the forward hadron matrix elements of the operators in Eqs. CEJ-© through 

^Es(P^\ QCD 0^%Jp,S) = 2(x n )u± d P, ... P , n , 

(p,S\^ CB d^lJp,S) = 2(x n ) Au±Ad S {tloP , 1 . . . Pfin} , (5) 

(p,S| QOD 0£<?.&>,S) = 1 ^(x n ) Su±5d S [ ^ l] ...p lin+l} , 

where p is the momentum of the hadron, Mh is its mass and S is its spin. The plus or minus signs in Eq. (J5J 
correspond to choosing isospin index or 3 respectively. The corresponding off-forward matrix elements are similarly 
related to Bjorken-x moments of generalised parton distributions (GPDs) which shall be discussed briefly below (see 
|10| for a comprehensive review). 

The hadronic matrix elements of the twist-two operators in Eq. © can be calculated using standard lattice tech- 
niques. Although a parametric form must be assumed in order to invert 0,0] the Mellin transforms, Eqs. such 
calculations will then lead to information about the parton distributions directly from QCD 1 . However, all lattice 
calculations are necessarily performed on finite volumes and at finite lattice spacings. Additionally, with current com- 
putational resources, statistically meaningful simulations can only be performed at quark masses, m q , considerably 
larger than those found in nature. These three restrictions have significant effects on calculations of twist-two matrix 
elements which must be taken into account if realistic predictions are to be made. 

Conveniently, the low energy QCD dynamics that these matrix elements characterize can be described using effective 
field theory. Standard chiral perturbation theory (%PT) as formulated in the infinite volume continuum allows 
systematic exploration of the quark mass dependence of low energy hadronic observables in the region where \p\ < 
A x where p is a typical momentum and A x ~ 1 GeV is the chiral symmetry breaking scale. Extensions to include 
finite volume (FV) and finite lattice spacing effects are also well developed (see Refs. ^jj and 0] respectively for 
recent reviews) as are the modifications necessary to treat valence and sea quark masses independently - quenched 
and partially-quenched xPT (QxPT and PQxPT) Q . In our study, we shall ignore the effects of the 
discretisation of space-time 2 (whereby our results will only be strictly applicable to lattice calculations in which a 
continuum extrapolation has been performed) and consider continuum partially-quenched chiral perturbation theory 
in a finite spatial volume of dimension L 3 . If the size of the box is large compared to the inverse pion mass (the 
lightest asymptotic state), M^L 3> 1, the power counting of infinite volume xPT (p-counting) applies and the 
necess ary modifications are easily made, replacing momentum integrals by sums over allowed momentum modes (see 
Refs. H3, [13, H2, E3, [S3, |iH HE E3, HE Ha H3, E3, for recent examples). On the other hand if M„L - 1, one 
needs to treat pion zero modes (components of the pion field with zero momentum) carefully since they correspond 
to vacuum fluctuations of order unity. In such a regime, modified power countings are required [33l l34l l35|. In this 
analysis, we will restrict ourselves to the region M^L ^> 1. 

Using the low energy effective theory, it is possible to compute the quark mass and volume dependence of hadronic 
observables such as the matrix elements of twist- two operators. For the most part, the quark mass dependence of the 



1 The reduced symmetry of the hyper-cubic lattice leads to lower dimensional operator mixing for twist-two operators with n > 3, and 
consequently calculations are only currently available for n = 1, 2, 3. 

2 The additional, Lorentz non-invariant contributions to unpolarised twist-two operators that must be included when the lattice spacing 
is non-zero have been considered in Ref. Illl . 
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various twist-two matrix elements has been studied extensively [36l Wh . l38l l39l I40L EH E3 ] and lattice data have been 
analysed assuming an infinite volume [36l Eiil Eij . However, the volume dependence of these observables has been 
ignored (with the exception of the matrix element of the n = helicity-dependent twist-two operator, the isovector 
axial coupling |30l l34jp in such analyses. Nonetheless, finite volume effects have been found to be important in 
many observables; here we investigate the effect they have on nucleon, and other octet baryon matrix elements of 
twist-two operators. 

In Section [n] we introduce aspects of heavy baryon chiral perturbation theory relevant for the analysis of twist-two 
matrix elements and define our notation. In Section lTTTl we discuss the twist-two operators in QCD and their matching 
in the low-energy effective theory and present examples of results for the quark- mass dependence of the nucleon matrix 
elements using two degenerate flavours of quarks. Full results in the two flavour partially-quenched case and results 
including the strange quark are relegated to Appendices iBl and lUl respectively. In Appendix iDl the quenched theory 
is discussed. In Section IIVI we discuss the general form of finite volume corrections to these matrix elements and 
make comparisons with available data. Section IIV CI discusses the complications that arise when non- forward matrix 
elements are considered and Section presents our conclusions. 



II. HEAVY BARYON CHIRAL PERTURBATION THEORY 



Heavy baryon chiral perturbation theory (HB^PT) was first constructed in Refs. [4rL Efil EH Esj . In current lattice 
calculations, valence and sea quarks are often treated differently, with sea quarks either absent (quenched QCD) or 
having different masses than the valence quarks (partially-quenched QCD) 3 . The extensions of HB%PT to quenched 
HB^PT [43 and partially quenched HB^PT |4ll |42| to accommodate these modifications are also well established 
and have been used to calculate many baryon properties |IllMIIllIlEilIHEl|53|. In this and the next sections, 
we will primarily focus on the two flavour partially-quenched theory; here we briefly introduce the relevant details 
following the conventions set out in Ref. [41|. We leave the three flavour and quenched cases to Appendices ICl and IDl 



Pseudo-Goldstone mesons 



We consider a partially-quenched theory of valence (u, d), sea (J, I) and ghost (u, d) quarks with masses corre- 
sponding to the matrix 



toq = diag (m u ,rnd,rnj,mi,mz,,mj) . 



(6) 



where 



m u d such that the QCD path-integral determinants corresponding to the valence and ghost sectors 



exactly cancel. 

The corresponding low-energy meson dynamics are described by the SU(4|2) PQxPT Lagrangian. At leading order 



f 2 f 2 
U = V str [Wrtdp-E] + A Vstr [m Q Et + m Q E] + a*0"$ o 0/,*o - rog$jj , 
8 4 

where the pseudo-Goldstone mesons are embedded non-linearly in 

'2 i $ 



with the matrix $ given by 



where 



$ = 



M 



( Vu 


7T+ 


J° 


L+ \ 




>]<: 


J- 


L° 


J° 


J+ 


V.j 


Y + 




L° 


Y ii 


m j 



M 



exp 



M x f 
X M 



Vu ^ 

k~ fjd 



f 



Xriu Xir+ X.JO XL+ 
Xtt- Xvd Xj- Xlo 



(7) 



(8) 



(9) 



(10) 



3 At finite lattice spacing, different actions can even be used for the different quark sectors (e.g., staggered sea quarks and domain wall 
(DW) valence quarks), but we do not consider this complication here. 
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The upper left 2x2 block of M corresponds to the usual valence-valence mesons, the lower right to sea-sea mesons 
and the remaining entries of M to valence-sea mesons. Mesons in M are composed of ghost quarks and anti-quarks 
and those in \ of ghost-valence or ghost-sea quark-anti-quark pairs. Because of the graded symmetry of the partially- 
quenched theory, the mesons in \ are fermionic. In terms of the quark masses, the tree-level meson mases are given 
by 



A 



(m 



Q) 



K) 



33 



(11) 



where Q = (u, d, j, I, u, d). 

The singlet field $o = str (<f>) / \[2 has mass mo at tree level. The terms proportional to a$ and mo in Eq. J7J) are 
only relevant in the quenched theory (see Appendix lD")l : in PQxPT and xPT the singlet mesons acquire large masses 
and can be integrated out. Furthermore, a* is suppressed by 1/N C and we set it to zero throughout. 



B. Baryons 

In SU(4|2) HBxPT, the physical nucleons (those composed of three valence quarks) enter as part of a 70-dimensional 
representation. This is described by a three index flavour-tensor, B 0,E3J3- The embedding of the physical nucleon 
fields into B and the symmetry properties of B are described in Ref. |41|. The A-isobar must also be included in 
the theory since the mass-splitting, A, between the nucleon and A-isobar is ~ 300 MeV, comparable to the physical 
pion mass (and less than pion masses used in current lattice simulations). The parameter A is assumed to be small 
compared to the chiral symmetry breaking scale. These fields are represented in a three index flavour-tensor T M (a 
Rarita-Schwinger field) transforming in the 44-dimensional representation of SU(4|2). 

The relevant part of leading-order Lagrangian describing these baryons and their interactions with Goldstone mesons 



L B = i (Bv ■ VB) — i (t% • VT^j + A 

+2a (BS^BA^) + 2(3 (BS^A^B) + 2H (t^S^A^ 



T A V B) + (BA V T U ) 



(12) 



where is the baryon velocity, is the covariant spin- vector |45j, |4j] and is the usual covariant derivative 

V»B = d"B + [V,B] and £>% = <9% + [V u , %] . (13) 
The vector and axial-vector currents appearing in the above expressions are given by 



v = - ( 



A» = - ( £0fft - gd^i 



(14) 



where ^ is defined in Eq. JHJ. The various Lorentz and flavour contractions (indicated by the parentheses) are defined 
in Ref. |4l|. In order that correctly describes the spin-3/2 sector, v ■ T = S ■ T = 0. 

In what follows, we will substitute a = ^gA + |<?i, P = §<7i — \9A, C — —gAN and H = <?aa since these correspond 
to the usual %PT couplings when the QCD limit, where mj — m u and m; = m<2, of the theory is taken. 



III. TWIST-TWO OPERATORS AND MATRIX ELEMENTS 



A. Twist-two operators in (PQ)xPT 



In the low energy effective theory, the twist-two quark bilinear operators in Eqs. 0-© match onto hadronic 
analogues constructed to obey the same symmetry transformation properties. In two flavour QCD, the unpolarised 
and helicity operators transform as either (3,1) © (1,3) (isovector) or (1,1) (isoscalar) of SU(2)lxSU(2)^. When 
one considers SU(4|2)£xSU(4|2)fl partially quenched QCD, there is more than one way to extend these operators 
I41II42I IH^j . Imposing super-tracelessness and the correct QCD limit in the valence sector, the most general extension 
of T3 (isovector) to the adjoint representation of SU(4|2)i ^ is 



t 3 = diag (1, -1, qj,qi,qk,qj + qi - <jfc) 



(15) 
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The freedom in choosing the values of the qi's can be advantageous in lattice simulations; certain choices of the q^s 
eliminate disconnected contributions (diagrams in which the operator is on a quark line connected to the external 
states only through gluons which are notoriously hard to compute |Hfij ) even away from the isospin limit. The 
non-uniqueness of the extension of Gell-Mann flavour matrices to PQQCD has additional consequences in that are 
discussed in Appendix IO 

For the isosinglet operator, the most convenient choice is 



f = diag(l,l,l,l,l,l) , 



(16) 



because it is purely in the singlet representation of SU(4|2). Any other choice, such as diag (1, 1, 0, 0, 1, 1) [which one 
might choose as disconnected diagrams would be absent], will contain contributions from other representations, and 
hence introduce additional low energy constants. 

The transversity operators in QCD are chiral-odd and belong to the representation (2, 2) © (2, 2). The most general 
choice for their extension to SU(4|2) PQQCD is 



t t = diag(l,yi,yj,y k ,yi,y m ) 



(17) 



For this operator disconnected contributions vanish as the matrix element involves helicity flip. Thus clean calculations 
of (x n )su and (x n )sd are possible. 

Based on these symmetry properties, at leading order in PQ^PT the hadronic operators that match onto those of 
PQQCD are 



0(A) 



it™): 



■n+1 p 



A« 4 



str 



£ + T A du . . . a Mn s + £f A d M0 . . . d Mn £t 



■ W/m-a ( 



T, 



traces , 



(18) 



0(A) 



Aa[ rA) %„ ■ ■ • v^_ t (BS^Bf^) + A/#*>» { *> . . . (BS^ n} f A + B) 
+ (1 - S A0 )Ac^v {ll0 . ..v, n _ ± [(T^b) + (Br| + T^ } ) 



traces , 



(19) 



na = Sa n v {tl0 . . . v M (BS a] &f£) + 8(3 n v { ^ . . . v M (BS a] f^B) (20) 



traces , 



where t a = | (^Ta£ ± ^ta^) an d t^* = i (^tt^ ± £tt£) , and the different Lorentz and flavour contractions 
(indicated by the parentheses) are given in Ref. ^l|. The super-script on the low energy constants (LECs; a^ rA \ 
Aa( rA ', etc.) in the unpolarised and helicity operators labels the chiral representation to which they belong; for 
A = 0, ro = s (singlet) otherwise r A = a (adjoint). In what follows, we take A to be either or 3. In QCD, the two 

different flavour contractions of the operators proportional to an A ^ and f3n A ^ (and their spin dependent analogues) 
are identical. 

There are additional classes of operators that formally enter these expressions at the same order but do not 
contribute to the next-to-leading order (NLO) matrix elements, i.e., their contributions to one-loop diagrams vanish; 
for example, 

w {,uo • • • *V„-2 ( r ^-i S ^} f A B ) » and v {^o--- v ^-i ( BS n n } T A B ) ■ ( 21 ) 
Such operators are omitted in Eqs. ^SJ-JSOJ. Also, NLO counter terms, such as 

(b[t a \m + )b) , where M+ = ^ (?m Q ^ + Zm Q £) , (22) 
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(a) 



(b) 



(c) 



o 



(d) 



(e) 



✓ \ 



o 



(f) 



(g) 



(h) 



O 

0) 



(j) 



FIG. 1: Diagrams contributing to nucleon matrix elements of the twist-two operators. The black square corresponds to an 
interaction from the strong Lagrangian and the gray circle represents an insertion of the twist-two operators in Eqs. 118H - 
I2UH . The thin, thick and dashed lines are 70-plet baryons, 44-plet baryons and mesons respectively. The first two diagrams 
represent the wave-function renormalisation and the remainder are operator renormalisations. Diagrams (e) and (f) are absent 
for the unpolarised, and isoscalar operator matrix elements, and diagrams in which the twist-two operator is inserted on a 
meson line are only present in the unpolarised case. 



are neglected in this work, since we are focusing on finite volume effects arising from one-loop diagrams at NLO. For 
the unpolarised isovector operators, these counterterms are explicitly displayed in Ref. |4l|; for the other operators, 
they are simple generalisations. 

Additional, higher-order operators arise when powers of the baryon velocity are replaced by derivatives, such 

as V{^ . . . ( ■fjd) fJiTl _ 1 ^BS^yT^ B^j. In the forward limit, these operators only appear in loop diagrams, so their 

contributions to matrix elements start at next-to-next-to-leading order (NNLO), therefore we do not include them in 
this work. 



B. Nucleon matrix elements 



The one-loop diagrams that contribute to nucleon twist-two matrix elements at NLO are shown in Fig. ^ The 
first two diagrams, (a) and (b), represent the wave- function renormalisation whilst the other diagrams are operator 
renormalisations. Diagrams (e) and (f) are absent for the unpolarised or isoscalar operator matrix elements as the 
transition between 70-plet and 44-plet baryon states changes spin and isospin. Finally, diagrams in which the 
twist-two operator is inserted on a meson line [(g), (h) and (j)] are only present in the spin-averaged cases. 

In Appendices FBI lUl andlDl we give the results for the independent matrix elements in SU(4|2) PQ^PT, for SU(6|3) 
PQxPT and for SU(2|2) quenched xPT in the isospin limit. As an example, here we present the SU(4|2) isospin limit 
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(m u — m^, rrij — mi), Qj — qi result for the nucleon matrix element of the isovector, unpolarised operators O^.J... 



(NP^JN) = -U p v Uo . ..vM2a^ - (3^) x (1 + (1 - 8 nQ )W su{m ) 



i(l - S n0 ) 
6/ 2 



8 2 I (a) °™ 
/!■ ■ //.., ^ g <l \ \ I 7n ~ 



' Van ■ ■ ■ V u 



3H(AU,A) +2H{M uj ,A) 



+a 



(a) 



8iX(M uj ) -3H(M uj ,0) 2 gi g A +3g 2 +3H(M V ,0) 2g 2 A + 2 gi g A + 3g 2 



+12 6 2 (gi+g A fH v ,{M„,0) 



4il(M uj ) - 12H{M uj ,0) g l9A + 3H(M„,0) [4 gi g A + g A ] 



+6 6 2 (g 1 +g A yH T1 ,(M 7V ,0) }, (23) 



where W' 



SU(4|2) 



is the nucleon wave-function renormalisation given in Eq. (|BlJ and 



S 2 = M 2 - M|. = A( 



(24) 



is proportional to the difference between sea and valence quark masses. The functions 1(M), H(M, A) and (M, A) 
are defined in Eqs. i|Al|l . i.YMi and i|34|) . below. Finally, Ub corresponds to the type B baryon spinor. To take the QCD 
limit, we would set 8 — > and j — > u. For eq uivalent choices of the f A our results reproduce those found previously 
for the unpolarised isovector operator [ill l42l| . One can also calculate the matrix elements of these operators in the 
A-isobar (and the N-A transition in the spin dependent cases). However since these are not stable particles in much 
of the region where X-PT converges, we do not present the expressions. 

In these results, the only effect of the diagrams in which the twist-two operator couples to a meson (diagrams (g), 
(h) and (j) in Fig.^l is to satisfy the number sum rule for the n = matrix elements, producing the <5„o factors in 
the above expression. For n > 0, these diagrams give sub-leading contributions, entering at 0(p n+2 ). The number 
sum-rule also fixes 



2a, 



(a) 


.,0) 



PS 



(a) _ 



0) 



3. 

3, 



t£° = 3, 



7o 



3, 



0. 



(25) 



and the n = low energy constants of the spin-dependent operators can be fixed in terms of the usual axial couplings 



2A. 



» 



A^ a) = 6g A , 



A 7 £ o) = 2g AA , 



A, 



» 



A< J +A/3« =2(g A + 9l ), A 7 £ 



(») 



= 0, 
0, 



AcW = - 



0. 



(26) 



IV. FINITE VOLUME CORRECTIONS 



A. General discussion 



In momentum space, the finite volume of a lattice simulation restricts the available momentum modes. Here we shall 
consider a hyper-cubic box of dimensions L 3 xT with T ^> L. Imposing periodic boundary conditions on mesonic fields 
leads to quantised momenta k = (ko,k), k — ^-j = ^-(ji, J2,jd,) with % € Z, but fco treated as continuous. On such 
a finite volume, spatial momentum integrals are replaced by sums over the available momentum modes. This leads to 
modifications of the infinite volume results presented in the previous section; the various functions arising from loop 
integrals are replaced by their FV counterparts. In a system where M^L ^> 1, finite volume effects are predominantly 
from Goldstone mesons propagating to large distances where they are sensitive to boundary conditions and can even 
"wrap around the world" . Since the lowest momentum mode of the Goldstone propagator is ~ exp(— M^L) in position 
space, finite volume effects will behave as a polynomial in 1/L times this exponential if no multi-particle thresholds 
are reached in the loop. 
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To investigate this behaviour, we consider the various finite volume sums occurring in the twist-two matrix elements. 
We first define 



£3 E I dk ° 



kuk v 



/2A 2 



(fc 2 - to 2 + ie)(k ■ v - A + it) 16tt 2 \ 3 



g ilv T{m, A) + v il v v Q{m, A) . 



(27) 



where the ultra-violet divergence has been subtracted in dimensional regularisation 4 with A = j 2 ^ — Je + log(47r) + 1 
(d is the number of dimensions). All finite volume sums that occur in the baryon wave function and operator 
renormalisations involving baryon propagators (diagrams (a)-(h) in Fig. ^ sunset-type diagrams) can be expressed 
in terms of Tim, A) and its derivatives. The tadpole diagrams in Fig. ^ are discussed in Appendix^] In the baryon 
rest frame where v = (1, 0, 0, 0), Poisson's summation formula allows us to decompose T into its infinite- volume limit 
and a volume-dependent part, 



T(m, A) = F(m, A) + F fv (to, A) . 
It is straightforward to show that the infinite volume piece is 



(28) 



F(m, A) 



16tt 2 





r , 2A 2 " 


A log 


/to 2 \ 


'10A 2 


4to 2 " 
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m 


b) + 






3 
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3 



A + l(A2_ m 2)3/2 lQg 

O 



A - VA 2 - m 2 + it 
A + VA 2 - to 2 + it 



(29) 



(fx is the renormalisation scale), and the finite volume corrections are given by 



i^ V (TO,A) 



E- 



12tt 2 ^ uL 



d\k\ 



mi»l im 



E 



-umL 



24ir ^ uL 



-A, 



where u = (u\, U2, 1*3) with m € Z, u = \u\ and 



A 



e (z2) [1 - Erf(z) 1 + (^—] 
\umL J 



9z 

T 



umL 



2 r 



39z llz 3 9z 5 
~64~ + ~~32 16" 



with 



\k\ sm(u\k\L) 
k\ 2 + m 2 + L 



A 



2z 2 ) e (z2) [l-Erf(z) ] 



z z 



e ( * 2) [l-Erf(z) ] 



O 



(urnLy 



(30) 



(31) 



umL 



(32) 



Higher order terms in the l/(umL) expansion in Eq. I|31|l are easily calculated. For convenience, we also define the 
functions 



H(m,A) 



dT(m, A) 

eK 



a r( a\ ^{m,A)-T(m,0) 
and K,(m, A) = 



(33) 



and their finite volume counterparts which we denote by the corresponding roman letter with the superscript FV as 
in Eq. e.g. K, -> K FY . 



4 It is important to note that because of the separation of scales, FV effects (infrared) are essentially independent of method chosen to 
regulate the divergent integrals (ultraviolet). Also, the results presented in this work are derived in Minko wski sp ace. We are free to 
work in Minkowski space since the sicknesses of quenched and partially quenched theories discussed in Refs. 20;, 5j| do not occur in our 
calculations. 



9 




■ A = MeV 

- A = 150 MeV 
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A = MeV 
A = 1 50 MeV 
A = 300 MeV 
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0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 

m (GeV) 



0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 
m (GeV) 



FIG. 2: Dependence of finite volume effects on the mass-splitting A in individual integrals/sums corresponding to diagrams 
(a) — (f) in Fig. The point m — 0.3 GeV corresponds to mL = 3.8. 




A = MeV 
A = 150 MeV 
A = 300 MeV 
A = 450 MeV 



0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 
m (GeV) 



FIG. 3: Dependence of finite volume effects in double-pole contributions on the mass-splitting A. Note that the scale here is 
ten times that in Fig. |5] The point m = 0.3 GeV corresponds to mL — 3.8. 



The function A represents the modification of the FV effects due to the mass-splitting A; in the limit A — > 0, 
A — ► 1. Figure shows the dependence of FV effects on the scale A for the functions if FV (m, A) and K FW (m, A) 
that arise from diagrams (a) — (f) in Fig. ^ It is clear that finite volume effects in individual diagrams involving 
a 44-plet are suppressed relative to those involving only meson and 70-plet baryon propagators, t houg h this can 
be compensated for by large coefficients. A very similar result was found in the heavy meson sector |23j where the 
contributions involving B* mesons are suppressed compared to those involving the B meson by the mass difference 
Ab = tub- — tub- However, the origin and behaviour of the mass difference in the current context is distinct. In 
contrast to the heavy meson case where As ~ I/thb arises from the breaking of heavy-quark spin symmetry and 
vanishes in the heavy quark limit, the mass difference A is generated by strong- interaction dynamics and remains 
finite in the chiral limit. Empirically, A is almost constant over the range of quark masses considered here. 

When one considers quenched or partially quenched theories rather than standard x^T, one expects somewhat 
larger finite volume effects because of the enhanced long-distance behaviour of double-pole structures in the singlet 
meson propagators of these theories |20j |. These double pole contributions are given by terms proportional to the 
functions 

H„'(m,A) = — — 2 — , /C n /(m,A)= — — - 2 — , (34) 

and the double-pole tadpole function (m) given in Appendix El [and their finite volume analogues constructed as 
in Eq. 128|) ]. From Figures |3 and it is clear that i? F , v (m, A) is about an order of magnitude larger than H FV (m, A) 
in accordance with expectations. 

In considering the magnitude of finite volume effects, the standard chiral power counting can be misleading; the 
FV effects of a diagram of a given order in the power counting may be larger than those of lower orders. For some 
generic observable, one may consider two contributions, C\ and C2, that enter at different orders, mi < m^, in infinite 
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volume xPT. As discussed above, the dominant finite volume effects in these contributions will typically be of the 
form SCi ~ (MffL) * exp(—M n L) when M n L ^> 1 and no multi-particle on-shell intermediate states can contribute. 
In some situations, the presence of additional meson propagators or other infrared enhancement in the higher order 
contribution (C2) can amplify its finite volume shift relative to that of the lower order contribution (Ci). For some 
(contemporarily relevant) choices of masses and volumes, the quantity 



{MM ,-H i'^'i >1; ,,,, 



??i2 — mi 

and the formally higher-order contribution will provide the dominant finite volume effect. In the current calculation, 
diagram (g) in Fig. ^ m which the twist-two operator is attached to mesonic propagators, may indeed fall into such 
a category. The finite volume corrections to these diagrams will be given by 

'4"»~«:,->^^<M.,0), (36) 
compared with those of the corresponding baryon operator diagrams (diagram (c) in Fig. ^) 

(a) 
f 



5lP~^H F v {Mv ,0). (37) 



From this we see that the ratio 



Sli n) (M v \ n 1 / M n L 



sip VA x y 4^ 



(38) 



where there is an undetermined coefficient involving ai?\ o4 a ' and other numerical factors that we assume to be 0{1). 
Whilst formally the magnitude of this ratio is indeed greater than unity for any n in the limit M n L — ► 00, both 5lP 
and 8 iff are exponentially suppressed. For realistic pion masses and volumes used in current lattice simulations this 
ratio is consistently smaller than one and the FV effects of diagrams (g), (h) and (j) in Fig. ^ can be neglected. The 
only exception to this is in the n = unpolarised matrix elements, the isoscalar and isovector quark numbers. Here, 
the volume dependence of diagrams in Fig. ^with mesonic operators exactly cancels that of those involving baryonic 
operators and wave-function renormalisations to give an overall result that is independent of the volume. 



B. Relevance to lattice data 



Lattice calculations of twist-two matrix elements have a long history, with the first calculations occurring in the 
1980s 57] . Over the last decade, a considerable effort has been made to study them in detail with major contributions 
from the QCDSF, LHP, RBCK and ZeRo collaborations. In Tabled recent simulation parameters are shown. State- 
of-the-art lattice simulations are beginning to enter the region of quark masses and lattice volumes in which the use 
of NLO chiral perturbation theory is justified (naively, this requires m 7T /A x < 1/3 and m^L ^4). At the moment 
however, there is little data in this region and a realistic fit of the combined mass and volume dependence in our 
PQxPT formulae (and thereby a determination of the LECs) is not possible; only general trends can be extracted. 

In order to address the expected size of finite volume corrections arising from our calculation, we first define the 
finite volume correction 

imam™ ( N \°\ N )T c - loop ~ (mo\N)T^ oop m] 

for each of the operator matrix elements we calculate. These corrections depend on a number of low-energy constants 
and couplings, some of which involve the A resonance. In principle, all of these parameters can be extracted from fits 
of the xPT forms to lattice data on nucleon matrix elements (thereby bypassing issues of the structure of unstable 
particles and transition matrix elements), however this is not practical at the present stage. Therefore, to fix the 

twist-two low energy constants (a„ , Sj n etc) we assume that large- N c relations [3^ amongst the parton distributions 
in the nucleon, A-isobar and A'-A transition are valid, leading to 

7 W = 2ai a) - pi a) , A 7 i a > = 1 (2Aa« - A/^>) , A<>> = \ ( 2 AaW - A/jW) , (40) 
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Group 


m T [MeV] 


Volume [fm 3 ] 


Notes 




Dynamical simulations (Nf = 2) 


QCDSF [58] 


560-940 


l.l 3 , 1.5 3 , 2.2 3 


Clover, a-0.08-0.12 fm 




LHP [59] 


340 


3.5 3 


Staggered sea, DW valence, a^ 


-0.13 fm 


LHP [5J 


340-730 


2.6 3 


Staggered sea, DW valence, a<- 


-0.13 fm 


LHP-SESAM [60J 


730-900 


1.6 3 


Wilson, a-~0.1 fm 




LHP-SCRI [60] 


480-670 


1.5 3 


Wilson, a~0.1 fm 




RBCK [£1] 


560-700 


1.9 3 


DW, a~0.13 fm 




Quenched simulations 


QCDSF [62] 


580-1200 


1.6 3 


Clover, a-0.05, 0.07, 0.09 fm 




QCDSF [63] 


310-1000 


1.5 3 


Wilson, a~0.09 fm 




QCDSF [64] 


440-950 


1.5 3 , 2.3 3 


Overlap, a~0.095 fm 




LHP [60] 


580-820 


1.6 3 


Wilson, a~0.1 fm 




RBCK [65] 


390-850 


1.2 3 ,1.6 3 ,2.4 3 


DW, a~0.15 fm 




ZeRo [66] 


750-910 


1.1 3 ,1.5 3 ,2.2 3 ,3.0 3 


Clover, a~0.093 fm 





TABLE I: Summary of recent lattice calculations of nucleon twist two matrix elements. Not all calculations involve the full set 
of twist-two operators. 



7« =a^ s) +/3i s) , A 7 i s) = Aa« + A/#> , Sa n = -A8p n = % n = ~5c„ . 

5 9 

The remaining LECs are not constrained by large- N c relations in QCD, and for want of accurate lattice data with 
which to fit them, we set fin'^ = c\4 , A/3^ a ' s ^ = Aa„' s ' and a^'^ = Aa^'^ — Sa n = 0. Throughout, we use 
/ = 0.132 GeV, and keep A = 0.3 GeV fixed independent of the quark mass. For the parameters appearing in the 
the flavour matrices T3 and fx, we set qj = q\ = qu = and yj — y^ — yi = y m = 0, and set yi to be either ±1. 
As discussed above, if one is using lattice data to determine the LECs, the g's and y's are fixed by the details of the 
lattice calculation. After making all of the above substitutions, the isospin limit results become proportional to the 
corresponding bare matrix elements and the finite volume effects given by Eq. I|39|l are easily studied. 

The axial couplings gA, gi, and g&N occurring in our results are the chiral limit couplings and there is some 
uncertainty in their values. We will fix gA = 1.3 (though even the chiral limit value of this is not well known [67|), 
|<?Aiv| = 1-5 and vary g± = ±gA- In the QCD limit, results are independent of g± since in this case it only involves 
couplings to the 77' meson which remains massive in the chiral limit and can be integrated out. 

Using these parameters, Figures 0] and [S] illustrate the typical size of finite volume effects in the various matrix 
elements. In Fig. QJ we consider a (2.5 fm) 3 box with a sea quark mass set such that the corresponding sea-sea 
Goldstone boson has a mass Mjj = 0.35 GeV and take g\ — gA- Fig. [5] is similar except here we take 171 = —gA to 
show the effect of this undetermined parameter. Variation of gA and <7a./v within reasonable bounds leads to similar 
modifications as those for varying g\. At the smallest pion mass in these plots, M^L ~ 3 and one must start to worry 
that infinite volume p-counting is no longer appropriate; at the largest pion mass in the plots, M 7r /A x ~ 1/2 and one 
must worry that the convergence of the chiral expansion becomes questionable. From these figures, it is nevertheless 
apparent that NLO PQxPT predicts finite volume effects in twist-two matrix elements that are generically < 5-10% 
for the range of masses and volumes studied here. However, there is some evidence that finite volume effects from 
higher orders of the standard chiral power-counting can be significant [68l l69| . 

Recently, staggered-sea, domain-wall-valence results have become available from the LHP collaboration |5jJ for 
very large volume (L = 3.5 fm) calculations of gA at a pion mass of 337 MeV. Also available are results on a 
somewhat smaller lattice (L — 2.6 fm) at the same pion mass. Although these two data points are consistent within 
their statistical errors (which will be reduced by ongoing calculations), their central values differ by ~ 15%. If we 
ignore the issues of non-locality due to the "fourth-root" trick used in calculating the staggered quark configurations 
and possible unitarity violations arising from the different valence and sea quark actions (which must vanish in the 
a — > limit that we have assumed), one can ask whether the NLO PQxPT formulae presented here can describe this 
dependence. To address this question, we consider the isospin symmetric QCD limit (m u = = rrij = mi) and 
define 

SgA = (A^ 3) iAor~ loop - (N\d ( ^\N)° L z^ oop (41) 



(A^ 3) |7V) 



tree — level 
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FIG. 4: Indicative finite volume effects in SU(4|2) matrix elements. The results in the first row are for the isovector unpolarised 
(left) and helicity (right) operators and those in the second row are similarly for the isoscalar unpolarised (left) and helicity 
(right) operators. The third row corresponds to the transversity "isovector" yi = — 1 (left) and "isoscalar" yi — 1 (right) matrix 
elements. In each plot, the solid curve shows the total result, whilst the short-, medium- and long-dashed curves correspond to 
the individual FV effects arising from diagrams (c)-(f), diagrams (a) and (b), and diagram (i) in Fig. Q In all of these results, 
we have considered a (2.5 fm) 3 box and set gA = gi = 1.3 and |<?ajv] = 1-5. M v = 0.25 GeV corresponds to M V L = 3.2. 



For this case, the LECs can be expressed in terms of the axial couplings through Eq. (|26[) and the FV shift, SgA, 
depends only on the pion mass, the volume and the chiral limit couplings gA, gAN and gAA- In Fig. we show 
the FV shift in gA that NLO xPT predicts at the LHP pion mass, M„ = 337 MeV. To illustrate uncertainties in 
the results, we vary the different axial couplings. In the central fits (indicated by the curves in the plot), we set 
gA = 1.3, gAN = —1 and g&A = —3 whilst the shaded band corresponds to 5gA for 1.0 < gA < 1.5, < |<7ajv| < 2 
and gAA — —3. Whilst, a shift of 15% between L = 2.6 and 3.5 fm is not predicted, the FV effects are substantial. 
However, without accurate knowledge of the chiral limit couplings, even the sign of the finite volume correction to gA 
is not well determined. 

As discussed in the previous subsection, in quenched lattice calculations, FV effects will be enhanced because 
of the double-pole contributions to singlet meson propagators. In Fig. [7| we plot the volume dependence of the 
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FIG. 5: As in Fig.gjexcept with 
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-g A . M n = 0.25 GeV corresponds to M n L = 3.2. 



polarised, isovector twist-two matrix elements in SU(2|2) quenched xPT (t ne analytic forms of these results are 
presented in Appendix lD|l . In contrast to partially-quenched %PT, in the quenched theory the LECs occurring in 
the Lagrangian and the twist- two operators are unrelated to those in standard %PT (though we denote them by the 
same symbols for convenience). To be definite, we choose Too = 0.7 GeV and the quenched operator LECs to satisfy 
Aa4 a ' = A(3n^ = 5A7i a ^ = 2Ac^ (as in the partially-quenched case), set the quenched axial couplings to gA = 1.3, 
|5aat| = 1-5 and let gi and 7 vary between ±gA as indicated by the shaded region. The curves in the figure correspond 
to gi = gA/2 and 7 = 0. As expected, the volume dependence here is enhanced over that in the PQxPT and xPT 
cases. 



C. Off- forward matrix-elements 



The off-forward matrix elements (in which the incoming and outgoing hadrons carry different momenta) of the 
twist-two operators correspond to moments of generalised parton distributions. Very little is known from experiment 
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FIG. 6: Finite volume effects in QCD calculations of gA at M n = 337 MeV (as appropriate for the recent LHP simulations 
[HsJ). L = 2 fm corresponds to M^L = 3.4. The shaded region corresponds to varying 1.0 < <?a < 1-5, < <?ajv < 2. 
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FIG. 7: Finite volume effects in the isovector, helicity matrix elements in the proton in SU(2|2) quenched xPT- The shaded 
region corresponds to variation of the axial couplings g\ and 7 between ±<7a with gA = 1-3, <?ajv = 1-5 and assuming large N c 
relations for the operator LECs. The point = 0.25 GeV corresponds to M^L = 3.2. 



about GPDs though major programs are underway at HERMES, Jefferson Lab and COMPASS to investigate them. 
As such, moments of GPDs are important quantities to extract from lattice calculations and much progress is being 
made in this direction |70j . It is therefore important to investigate the quark mass dependence 5 and size of finite 
volume effects in these calculations. Here we shall only discuss the novel features that appear in regard to finite 
volume effects when off-forward matrix elements are considered. A full analysis of the low-energy behaviour of these 
matrix elements will be given elsewhere [73) . 

To be specific, we shall consider the proton matrix elements (p'| c ^ CD C'^.' ) .. At „ \p) in which four momentum = 
(p' — p)^ (with \q 2 \ <C A^.) is injected through the twist-two operator. The analysis of these matrix elements is 
significantly more complicated than in the forward limit. This is primarily because the number of possible independent 
tensor structures in the matrix element grows with n; for example, 



(p'I qcd <^L>> = u p {p') 



7 ^V 1 . . . g*p"*+* . . .p^U^V) (42) 



5 Refs. El El address this issue for the infinite volume n = 1 matrix element relevant for the spin content of the proton in SU(2) HBxPT. 
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.0.1 



2m 



.p^B^tf) 



U p (p) , 



where p = \{p' +p). For each of the coefficient functions A^(q 2 ), B)^.{ci 2 ) and C^(q 2 ), there is an independent 
finite volume xPT expansion. 

Another complication enters when one considers the operators that match onto the twist-two operators in the low 
energy effective theory. The presence of the new scale q means that considerably more operators must be included in 
Eqs. (|18fl - H2U() ; f° r example, the term proportional to a4 in Eq. (|18|) is replaced by 



£< 



,{Mo 



. «"* (-id)^+ x . . . {-id)^ - tr + a%±\ [(-idY 10 . . . (-id)^ - tr 



(43) 



Note that only terms with an even number of derivatives contribute here. Each LEC in the forward case is replaced 
by 0(n) LECs. Additionally, the vector q allows more tensor structures to enter and we must also include 



n 
.7=0 



(-id) 



i-idy 



,- lf » [ s ^h s - (- 

V M 



id)} 



,B)-ix 



(44) 



since S-q ^ 0. When we also take into account the 44-plet fields many more structures are possible since q-T ^ 0. 
Even for the n < 3 matrix elements that have been calculated on the lattice, a large number of LECs need to be 
determined. This makes a reliable extraction of the physical matrix elements from finite volume, unphysical mass, 
lattice calculations a challenging proposition. 

In terms of FV effects, the modifications for the off-forward case are relatively simple and it is worthwhile to 
investigate them in some detail. There are essentially two classes of diagrams: ones where the twist-two operator 
injects momentum into a (heavy) baryon field (e.g., (c) in Fig. and ones where a meson receives the additional 
momentum (e.g., (h) in Fig. The former class is relatively uninteresting since for heavy fields, derivatives in the 
twist-two operators pick out only the momentum transfer between the external states, q, and can therefore be factored 
out of the integral. For this type of diagram, finite volume effects will be the same as those in the forward limit as 
we are free to work in the Breit frame where q ■ v = 0. 

For diagrams in which the twist-two operator is on a meson line, the situation is different since the derivatives in 
the operator can result in powers of the integration momentum. The relevant integrals are of the form 



-y 



L 3 



dk S ■ {-k)S -(k + q) . . . (-fc)^ (k + q)^ ...(k + q)^ 



tr 



2tt 



(k ■ v - A + ie) (k 2 - m 2 + ie)((k + q) 2 
+ a)S-(£ + b) 



(£ 2 -M 2 ) 3 



m 2 + ie) 

(-l) j \{i + a)^ ...(£ + o)^ (£ + fc) w+1 ...{£ + b)^ - tr 



(45) 



after introducing Feynman and Schwinger parameters and shifting the momentum integration k — >£ = k + xq — Xv. 
Here a — —x q + Xv, b = (1 — x) q + Xv and 



M = M 2 (x, A, m%q\ A) = m - x(l - x)q 2 + A + 2AA . 



(46) 



The trace subtractions prevent any of the £ A ' i 's arising from the operator from contracting with one another, conse- 
quently the non- vanishing scalar integrals/sums whose finite volume effects we are interested in will be of the form 



dx 



(IX 



1 



(£ 2 Y 



L 3 ^ J 2ir (I 2 - M 2 ) 3 



(47) 



where r — 0, 1 or 2. 

Without going into further details of the tensor structure j?^ , it is already clear from Eq. (|46|l that the effect of 
the momentum injection on overall finite volume shifts is very similar to the effect of the N-A mass splitting. Since 
x(l — x) > 0, when space-like momentum (q 2 = — \q^ 2 < in Minkowski space) is injected, FV effects are suppressed as 
the meson receiving the momentum injection moves further away from its mass shell. However, if time-like momentum 
is injected the situation is more complicated. Provided the virtual particle cannot reach its mass shell, finite volume 
effects are enhanced over the forward case but will still remain formally exponentially suppressed. However, if the 
injected energy-momentum is enough to put the intermediate particles on shell, it leads to a cut in Minkowski space 
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in infinite volume. In this case, finite volume effects are only suppressed by powers of 1/L in QCD. In (partially) 
quenched QCD, volume corrections for isoscalar twist- two matrix elements may be proportional to positive powers 
of L whereby the infinite volume limit will be undefined. This suppression of finite volume effects with space-like 
momentum injection and enhancement in the time-like case (which is relevant for twist-two matrix elements between 
states of different masses, e.g. N — > A transitions) will occur in hadronic form-factors that are specific cases of 
twist-two matrix elements. 



V. CONCLUSIONS 



We have studied the matrix elements of twist-two operators that determine the moments of the unpolarised, helicity 
and transversity quark distributions to NLO in (partially) quenched chiral perturbation theory in both infinite and 
finite volumes. We have performed our calculations in Nf — 2 and Nj = 2 + 1 partially quenched heavy baryon chiral 
perturbation theory and also studied the SU(2|2) quenched theory. These results will be relevant for extrapolations 
of lattice calculations of these matrix elements in the proton and other octet baryons (e.g., the A hyperon |74|'). 

We have focused primarily on the effects of the finite volumes used in lattice calculations. Without accurate data 
in the chiral regime with which to fit the various low energy constants on which the results depend, it is difficult to 
be specific, however it is clear that for most current simulations FV effects are not negligible. For typical full- or 
partially-quenched- QCD calculations, they are < 5-10% but may be significantly larger in quenched simulations. 

In the case of the off- forward matrix elements relevant to generalised parton distributions, we have not presented 
full results for arbitrary moments |73j| . However, we have analysed the finite volume effects in these matrix elements. 
We find that they should decrease with respect to the forward matrix elements if space-like momentum is injected. 
On the other hand if time-like momentum is transferred, finite volume effects will be enhanced; in QCD they may 
become only 1/L suppressed, and in (partially) quenched QCD finite volume isoscalar matrix elements may even be 
proportional to powers of L. 
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APPENDIX A: TADPOLE INTEGRALS AND FINITE VOLUME SUMS 



The sums that appear in tadpole diagrams, after subtracting A, are 

x 1 f dhn i m 2 - , . . 

X m = 73 E / ^TT2 2-^ + 77T^ A > A1 

L- 3 J 2n k z — m + it 1oit z 



k 

and 



1 1 1 } _ dl(m) ^ 

(k 2 — m 2 + it) 2 167T 2 dm 2 

k 

Using Poisson's summation formula, it is straightforward to show that 

T(m) = I(m) + I FV (m), (A3) 

where 

T , n 4-d f d d k i m 2 - m 2 , fm 2 \ 

I{m)=^ d ——- A=-— log — , A4 

J (271-)° k z — m z + it I07H I67H \ M / 
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is the infinite- volume limit of I{m), and 



mL>l 



-y 



TO7T / 1 



47T 2 ^— ' V 2uZ, V ul, 



e -umL ■ j 
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8umL 128(umL) 2 



O 



(umL)' c 



(A5) 



APPENDIX B: RESULTS FOR SU(4 2) PQxPT 

In this section, we present the results for twist-two matrix elements in the isospin limit in SU(4\2) PQxPT. The 
various masses and the mass-splitting 5 are defined in Sections [H] and IIIII 
The nucleon wave-function renormalisation is 



Wsum = ^p{H(M*,0) {-5gj-4g igA + g 2 A )+H(M u] ,0) (5gf + Ag x g A + 8 , 

+4 (H(M W , A) + H(M u j, A)) g\ N + (-6 g[ -12g l9A -6 g\) S 2 H v/ (M 7V ,0)} 
The isovector, unpolarised nucleon matrix element is 

(Npfl.^N) = hjNV^ . ..v^U N (2a^ - ffl) x (1 + (1 - S n0 )W su(m ) 
12/2 — u nv^ . ..v^Un < -g AN I TA J 3- J x 



3 W(M„ A) (-4 + ^ + + H(M uj , A) (8 + 3 y + 3 <r) 



(Bl) 



- 4iX(M ttf ) (-4 + 3© + 3©) +4iX(M 7r ) + ©) 

+3 ( - H{M uj ,0) \Ag\ { qj + qi) + 2g igA (2 + q 3 + qi ) + g\ (6 + gj + g ; ) 
+W(M OT , 0) [4.g^ (1 + qj + qi) + 2 9l g A (2 + qj + qi) + g\ (6 + q 3 + qi 

+8S 2 (.g 1 +.g j4 ) 2 ^(M 7r ,0) 



4 i l(M uj ) (2 + 3 qj + 3 qi) - 4 i l(M n ) {q 3 + g; ) 

F3 (W(M UJI Q) ff i [-8 ffj4 + 3 pi (eft + gj)] + H(M„,0) [8 gi g A + 2g 2 A -3g 2 1 (q 3 + ©) 



+4<5 2 (.g 1 + « ?j4 ) 2 ^(M x ,0)) 



(B2) 



The isovector, helicity matrix element in the nucleon is 



16 /2 



12 /2 U NV{no--- v ^-i S Hn} U N <j y y I SAN A<4 Q) X 



/C(M W ,A) (8 5j4 + .gi (2-3g,-3g,)) +/C(M ui ,A) (8 5A +.gi (-2 + 3g, + 3g ; ) ) 



20 
27 



5ajv 



3 H(Jlf», A) (-4 + qj + q t ) - H(M uj ,A) (8 + 3 g, + 3 g,) 



(a)' 
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+Aa^ 



12 i 1(M W ) (q 3 + qi ) - Ul(M uj ) (-4 + 3^+3 qi ) 

+H(M uj ,0) (lg A { qj + qi ) + 2 9l g A (2 + q 3 + g,) + y 2 (6 + 9j + g,) ) 
—H(M n , 0) (4 5 i (1 + g^ + g,) + 2 fl i 5 a (2 + <& + ®) + g\ (6 + gj + ®) ) 



-8<S 2 ( 5l +.g A ) 2 H v (M T ,0) 



o) 



12 »Z(M W ) (gj + ®) - 4iZ(M uj ) (2 + 3 qj + 3 ®) 
+W(M uj , 0) 5i ( - 8 g A + 3 .gi fa + g ; ) ) 

+W(M„ 0) (8 9l g A + 2 - 3<tf ((& + ®) ) + 4 <5 2 (51 + . 9j4 ) 2 (M ff , 0) 

and the isoscalar, unpolarised matrix element is 

W&Ujtf) = ■ • -v^ n U N (a^ +/30l s )) x (1 + (1 - <5 Il0 ) W S[/(4 |2)) 



(B3) 



+ Jp U " v w ■ ■ ■ v^U N (l - S n0 ) 1 4 (W(M T) A) + H(M uj , A)) ,g| w ^ 
+ (^ s) +/3i s) ) 



( S ) _ ^_ 



W(M T ,0) (5 5 ? + 4 5 i 5 A-.gA) 



- H(M uj , 0) (5 .g 2 +4g l9A + 8 g\) + 6 5 2 ( 5l + .g A ) 2 (M„ , 0) 

The isoscalar, helicity matrix element is 

(N\0^JN) = U N v^ ...v fln _ 1 S flny U N (Aa^+A^) x (l+W S um) 
i ^ „ f 40 

"3 



(B4) 



Act. 



^ {w > ■ • ■ v^S^Us I y (W(M W , A) + H(M u j, A)j g 2 AN |^A 7 „ 
2H(M V ,0) (5 3 2 +4. 9i5 a- 3 2 1 ) 



(■•0 



12/2 
+ (Aa« + A/# 

+2H(M ttj , 0) (5g 2 + 4 3l 5A + 8 3 i) - 12 S 2 ( 9l + g A f H V ,(M V , 0) 

Finally, the transversity matrix elements are 
(N\Ol..., na \N) = 

\u N v^ a . . . v^S^Un^ + yi )8a n + 2(1 + 2y i )5(3 n ) x (1 + W su(m ) 



(B5) 



* T7 n „ 16 ^ 

+ T272 •■• v m s h*] Un \ Y\l3 9AN 



K{M uj , A) (4 3A (1 - yi ) + 9l (-4 - 2 y< + 3 y, + 3 y fe ) ) 



+/C(M 7r , A) (4.g A (1 - Vi ) + 9l (4 + 2 y 4 - 3 Vl - 3 y m ) ) 



Sc„ 



20 
27 



9 AN 



H(M uj ,A) (10 + 2 y, + 3 % + 3 y fc ) + 3 H{M n , A) (6 + 2 y, - y, - y m ) «5 7 „ - 



3 <5(7 r 



+(5a„ 



2 2 [Z(M u j) (10 + 2 y, + 6 % + 6 y k ) + Z(M X ) (1 + 5 y, - 6 Vl - 6 y m )] - 4 (5 + Vl ) S 2 Z,, (M x ) 
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+5f3 n 



+H(M uj , 0) (gl (7 + Vl + y 3 + y k ) + 2 9l g A (2 + y 3 + y k ) + Ag\ (1 + w + y 3 + y fc )) 
+W(M w ,0) (<?i (3 + 7 Vi - iyi - Ay m ) -2 9l g A (2 + Vl + y TO ) - g{ (7+ y< + j/j + y m )) 

-2 ( 5l + g A f (5 + W ) <5 2 Hrf{M v , 0) 

4 i [X(M U , ) (2 + 4yj + 3jy + 3y k ) + J(M W ) (2 + ^-3^-3 y m )] - 8 (1 + 2 W ) <S 2 I,, (M T ) 

+H(M uj , 0) (8 si g A + 8^(1 + w ) + g\ (2 + 2 w + 3 Vj + 3y k )) 
+H(M„ , 0) (2 g\ - 8 <?i <m y< - g\ (2 + 2 w + 3 y x + 3 y m )) 

-4 ( 5l + 5A ) 2 (1 + 2 Vl ) S 2 H v , {M n , 0) 



(B6) 



APPENDIX C: RESULTS FOR SU(6|3) PQxPT WITH m u = m d ± m s 

In three flavour QCD, one seeks to determine the up, down and strange quark, unpolarised, helicity, and transversity 
distributions in the octet baryons. Consequently, the goal of lattice calculations is to determine the corresponding 
twist-two operator matrix elements for each of these flavours in the octet baryons. The SU(6|3) results presented in 
this appendix will be relevant for chiral and infinite volume extrapolations of lattice calculations of moments of the 
the strange-quark distributions in the nucleon 6 and the various parton distributions in, for example, the A hyperon 

If we consider the extensions of three-flavour QCD to partially quenched theories, we are naturally led to SU(6|3) 
HBxPT. The Lagrangian of this theory is very similar to that described in Sec. ^ with some simple modifications. 
Obviously the meson field $ is enlarged, becoming a 9x9 matrix encoding the 80-plet of pseudo-Goldstone mesons. 
The octet baryons are now embedded in a 240 representation and the decuplet baryons in a 138 representation. 
Additionally, the couplings a, (3 and C in Eq. (jT2"|) are replaced by a — > |Z3 + 2F, f3 — > — f-D + F and C — > C 
so that the nomenclature is the same as in SU(3) xPT. For definiteness, the quark masses we consider are tuq = 
diag(m M , mj, m s , nij, mi, rrir, m u , md, m s ) and Q after Eq. I|ll|l is replaced by Q = (u, d, s, j, I, r, u, d, s). Further 
details are given in Ref. |42fl . 

To calculate the independent moments of the PDFs in three flavour QCD, one constructs three independent flavour 
combinations of operators. The standard choice is 



• -£V}{ 1 ' A 3,A 8 }V'- 



(CI) 



where T = 7, 775, u represents the appropriate Dirac structure and the Xi are the usual Gell-Mann basis for 
SU(3). The unpolarised and helicity operators are in the singlet (1,1) or adjoint (8,1)©(1,8) representations of 
SU(3)ixSU(3)fl. From a lattice practitioner's point of view, A3 is somewhat special since in the limit m u = m^ there 
are no disconnected contributions to matrix elements of such operators. On the other hand, both the singlet and As 
operator require such contributions and no choice of flavour basis can ameliorate the situation. In partially quenched 
QCD, there is again freedom in the extension of the above QCD operators; a natural choice with a smooth QCD limit 
is 



A 3 = diag(l, -1, 0, 1,-1, 0, 1, -1, 0) and A 8 = diag(l, 1, -2, 1, 1, -2, 1, 1, -2) . 
In our results, only a single adjoint representation operator is presented, corresponding to 

X Adj = diag(l,gi, -1 - q x , l,q 2 , -1 - q 2 , l,£?i, -1 - qi) , 



(C2) 



(C3) 



and to determine A3 and Ag, we set q 2 = q\ = Tl- Keeping q 2 qi, will allow SU(3) breaking effects to be analysed 
in detail. However, the singlet operator is uniquely defined 



A = diag(l,l,l,l,l,l,l,l,l) 



(C4) 



One can also study strangeness in the nucleon using two-flavour xPT J7.^ thereby sidestepping issues of the slow(er) convergence of the 
chiral expansion around the physical strange-quark mass. 
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and disconnected loops involving sea quarks are unavoidable. 

The transversity operators in QCD belong to the (3,3) (3,3) representation of SU(3)i xSU(3)# irrespective of 
the choice of flavour structure in Eq. (|C1J| . For the partially quenched QCD extension of these operators, we choose 
the operators built from 



A T = diag(l,yi,y 2 , 0,0,0, l,yi,y 2 ) 



(C5) 



and then setting {yi, y 2 } = {1, 1}, {—1, 0}, {1, —2} gives the required flavour combinations. 

In this section, we give results for the matrix elements in the proton, A , S + , S~ as well as the S°-A° transition. 
Other octet matrix elements are simply related to these by isospin symmetry In the results of this section, the 
low-energy coefficients (al^\ A-jn \ etc) occurring in the SU(6|3) versions of Eqs. H18[) - l|20() are different from those 
in SU(4|2). With this caution, we use the same notation. 

The SU(6|3) tree level meson masses, M ss , M sr , M ur and M S j, are defined through Eq. (|TT)l . 5 is defined in Eq. 



S 2 = M 2 S - M 2 r = X(m s - m r ) 



and 



M z x = - 



Ml + 2Mt. - 2 [5 2 + 25 



(C6) 
(C7) 



The QCD limit is easily recovered, taking S — > 0, S — > 0, j ' — > u, I — > d and r — > s. To make the presentation succinct, 
we define the following ratios, 



21 : 



S 2 (m* - Ml + 26~'< 



Ml - Ml + 5 2 + 2S 2 



(C8) 



23 



{mi 



— Mr 



25' 



25 A 



<t = 



2(Ml - Ml + S 2 + 2(5 2 ) 2 

J 2 (M 2 -M 2 s -26 2 ) 
l M2-Ml-25 2 -45 2 ' 



(C9) 



(C10) 



M 2 - Mi 



2S 2 Y 



Ml - Ml - 2S 2 - AS 2 



(Cll) 



s 2 


[m 2 -M 2 s + 25~ 2 ) 




{Ml - Ml) 


f K Ml - Ml + S 2 + 2~8 2 ^ 



(C12) 



and the functions. 



3 = 2z- 



~S 2 {Ml -Ml -25*) 



(Ml -Ml) (Ml -Ml- 2P -4* 



Quu — 2G?)',tm + ^luu + (1 — 25)2X) 

Q ss = £X n ,, ss +m ss + (l-®)lx, 



(C13) 



(C14) 



Run = a?V,» u + 23W UU + (1 - 23)W X , 
R ss = <mr,', ss + ®H ss + (l-®)H x , 



21 




22 



+2H ujr (D 2 + 3 F 2 ) - 2H UU (D 2 - 3 F 2 ) - 2H US (D 2 -6DF + 3F 2 ) 



-9i 



R ss (D - F) 2 + 4 R us F(F - D) + 4 R uu F 2 



(C24) 



for the S baryons, and 

;(H) 



S£/(6|3) 3 y 2 



5 H us + H sjr + H ss + 2 H ujr - 2 i (T ss -2T US + T uu ) 



C 2 



sjr ' "ss 1 "ujr 

-3 3H ujr {D - F) 2 + 2H S3r (D 2 + 3F 2 ) - 2H US (D Z - 6DF + 3F 2 ) - 2H SS (D 2 - 3F 2 ) 



-9i 4 R ss F 2 — 4 R us F(D - F) + R uu (D — F) 2 



(C25) 



for the S baryons. 



2. Isovector unpolarised matrix elements 



(Pl^o-^nb) = g X ^ 1 + W SC/(6|3)( 1 - (5 ™o)J 

(a)' 



+ 



°n0) 1 



f 2 



3We„+W^ r ) (5 + gi) + 3(l + (fc)« 



(C26) 



C 2 



i(-5 - gi ) Z^ r + (1 + g 2 )J ujr - i H ujr (^3D 2 -6DF + 5F 2 + (D 2 + F 2 ) qi 
+ AiHuu D{D- 2F - F Ql ) - i(l + q 2 )(D 2 + 3F 2 )H ujr - ^R UU (D - 3F) 2 (5 + Ql ) 



— (1 + 2«i) l ujr + (1 + <&)2^> - 2 i W„,> ( F 2 + F 2 (l + 2^) 



+ 4iH uu D(D-2F- F qi ) - 3i(l + g 2 )(r> - i^Wujv - - 3F) 2 (1 + 2 9l ) 



(A°I<...,JA°) - - U ^---Y A0{1 + qi) («M-2#>) x (l + <j (6 | 3) (l-M) 



+ 



[7 A ou Mo . . .^ n [/ A o(l + gi)(l - (5 n0 ) 



/ 2 



(a)' 



c 2 



(a) 



+ ^W uj r (7D - 9F) (£> - 3F) - ^H sjr (D + 3Ff 



■zH,,,. (ID 



12DF + 9F 2 

i 



: H US (5D-3F)(D + 3F) 



-q(D + 3F) 2 H S3r - ^q(5D 2 - 6DF + 9F 2 )li ujr 





+ 4 ( R ss (D + 3F)' Z - AR US (2D - 3F)(D + 3F) + AR UU (2D - 3Ff 



12 

- 1 ( 2X uj + l ur ) + q± ujr + l -U ujr (ID 2 + 6DF - 9F 2 ) - ^H sjr (D + 3F) 2 



-^H uu (7D 2 -12DF + 9F 2 ) 



-H US (13D 2 - Y2DF + 9F 2 ) - 3iq(D - F) 2 TL n 
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R ss (D + 3F) 2 - 4R US (2D - 3F)(D + 3F) + 4R UU (2D - 3F) ; 



(C27) 



U^ + v^ . . . iv„[/ s+ ( (-4 + qi )a^ + 2(1 + 2 qi )(3, 



x (i+w^a-M) 



+ 



U-z+v^o ■ ■■VnJJ-z+il - 5„o) 



f 2 {9\J n 3 

_ + 2 + Hujr] Q1 + O- + Q2) (ftujr + 2 ^>) 

-2i (r ss -2T us +T ulI )(l- gi )V 2 



(a) _ 



(a) 



- - 5 T<y> + I sjr (l + qi) + (1 + q 2 ) 



" 2 '^ u i r ' 
5i 



(-D -F)(D + 3F) + 2(D 2 + F 2 ) qi 

2 iH uu (AF 2 + (D-F)(D + F)q 1 



-H sjr (D-F) 



+iH us 4(£> 2 ~DF + F ) + (D z + ADF - F z ) qi 



--(1 + <z 2 )(F 2 + 2DF + hF 2 )U ujr - -(1 + fc)(£> - F) 2 W sjr 



F ss (F - F) 2 - 4F US F (D - F) + AR UU F 2 ( Ql - 4) 



-iW ujr ((£> - F)(£> + 3F) - 4F 2 

-iW 5jr (D — F) 2 + 2i Hun (d 2 + F 2 + 2F 2 qi 

+2 i H us (d 2 - 2DF - F 2 + 2(D - F)F q^j 
-i(l + q 2 )(D - F) 2 

R ss (D - Ff - m us F{D-F)+ 4R UU F 2 I (1 + 2 Ql ) 



Us-v^.-.v^Us- ((5 + 4 qi )a ( n a) -2(-l + qi )l3 { n a) ) 

^ '- x ( 1 

6 V 



W SX/(6|3)( 1 - 5 no) y 

H% r + Hi (2 + qi ) + Ut (13 + 11 91) 



+ - 



Us-v^ . ..v^Us-il - 5 n0 ) 



P 



- - 7 (a) - 
9 7 " 



(C28) 



+4 U% r (1 + qi) - (1 + Q2 ) (2 + wA r ) - 2z (T ss - 2 T MS + T„„) (2 + 9l ) J C 
qi + 5I s j r (l + qi) + (1 + 92) 



+^n ujr (D - F) 2 (1 + qi) + ^W sjr (^F 2 - 2FF + 5F 2 — (D — F)(D + 3F)qi 
-m ss (d 2 - 5 F 2 - 4 F 2 ^ - ^(3D - 5F)(D -F) + A(D 2 -DF + F 2 ) q^j 
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--(1 + q 2 )(D - F) 2 H U j r - -(1 + q 2 )(D 2 + 2DF + 5F 2 )H S 



+ 1 ( 4R SS F 2 - 4R US F(D-F) + R uu (D - F) 2 ) (5 + 4 Ql ) 



- - 2T ujr qi + T sjr (l + qi) + (1 + q 2 ) 
i H sjr ( (D + Ff + [D - F) (D + 3F) qi J - 2 i H ss ( D 2 - F 2 + (D 2 + F 2 ) qi 



-2 i H us (d 2 - 4DF + F 2 + (D 2 - 2DF - F 2 ) q^j 

+ q 2 )(D - F) 2 (2H ujr +n sjr ) +iH ujr ({D - F) 2 (1 + gi ) 
+ (m ss F 2 - 4R US F(D-F) + R uu (D - F) 2 ) (1 - qi ) J 



(C29) 



<S°|0#...,JA°) 



+ 



t/go^o ■ ■ -v^Ua (-1 + gl) 
4V3 

U^oV^a ■ ■ .V^UAoj-l + qi)(l - S n0 ) f i I ( Q ) _ On 



(<#> - 2/?i a )) x [l + \ (w$ (6|3) + W$ (6|3) ) (1 - 5„o) 



V3/ 2 



- ; ( 2H^ U + H^s + 'H^,r I C 2 



+a 



(a) 



^J„ jr + ^Hujr (D - 3F)(3F -F)- l -H sjr (D - F)(D + 3F) 



-iH uu D(D + F) — iH us D(D - F) 

- l - (r ss (D -F)(D + 3F) - 4R US (D 2 - 2DF + 3F 2 ) - 4R UU F {-2D + 3F) 



(a) 



- ^ujr + (D - 3F)(D + F) + -H sjr (D - F)(D + 3F) 

-2tH uu D{D - F) — 2iH us DF 



+ ^ ( R ss (D -F)(D + 3F) - 4R US {D 2 - 2DF + 3F 2 ) - 4R UU F {-2D + 3F) 



(C30) 



3. Isovector helicity matrix elements 



(p\oit-, n \p) 



U P V{^ ■ ■ ■ v^S^yUp^ + q^Aa^ + 2(1 + 2gi)A/&°>) 



6 



(l + W^(6|3)) 



+ 



U p v {ll0 . ■ ■v^ n _ 1 S f , n} U p f hi 



+ ■ 



f 2 
Si 12 



Act, 



3 (5 + gi ) + 3(l + g2 )W 



9 V 3 



JC UU {D + 3F)(-1 + gi ) - /C„ jr (3F> — F — 2F qi ) + 3(1 + «fc)(D - F)/C u 



(C31) 
C 2 

CAci a) 



+Aa 



(a) 



1 



(5 + qi)l ujr + (1 + q 2 )I, 



II f I ~r ^ H*i ;l~ 



4i 

T 



3 F> 2 - 6 D F + 5 F 2 + (F> 2 + F 2 ) gi 
W„„ D (D — 2F — F Ql ) + 1(1 + q 2 ){D 2 + 3F 2 )H ujr + ^R UU (D - 3F) 2 (5 + gi ) 



+A/3( a ) 



1/ 



i(l + 2gi) J ujr + (1 + g2 )F ujr + y H u jr (d 2 + F 2 (l + 2 qi ) 

1 



- H uu D (D — 2F — Fq x ) + *(1 + g2 )(F. - F) 2 H ujr + ^Ruu(D - 3F) 2 (1 + 2 9l ) 



25 



<A«...,JA°> = _ ^ {w -^-^> yA ° (1 + gl) (a«M-2A#>) x (l + W$ 6|s) ) 
J7 A ov {Alo . ..w jUn _ 1 S' jUri} C/ A o(l + gi) 



+ - 



/ 2 



54 7 " 5 



Ai 2 

i 

4 



IC US (D - 3F) + 2/C„„ D + K ujr (D + F) + 2q(D - F)t 



CAci a) 



l ujr + 21 sjr + 2q(i ujr +l sjr ^ j --^H ujr (7D <)F)(D- IF) 



+^n sjr (D + 3F) 2 + -^-Huu {ID 2 - 12DF + 9F 2 ) + -^H US (5D - 3F)(D + 3F) 

oD lo lo 



+— q(5D 2 - 6DF + 9F 2 )H ujr + —q(D + 3F) 2 H sjr 

lo 



I 

T8 
1 

36 



R ss (D + 3F) 2 - AR US (2D - 3F)(D + 3F) + AR UU {2D - 3F) : 
- X -X U]T + ql ujr - ^H ujr (ID 2 + 6DF - 9F 2 ) + ^H sjr (D + 3F) 2 



~yH uu (ID 2 



12DF + 9F 2 



-Hus(l3D 2 - \2DF + 9F 2 ) + iq(D - F) 2 H V 



+ — ( R ss (D + 3F) 2 - 4R US (2D - 3F)(D + 3F) + 4R UU (2D - 3F) : 



(C32) 



t/ E+ ^ . . . v^S^yU^ ((-4 + qi )Aa { n a) + 2(1 + 2q 1 )A& ) ) 



6 



x 1+W 



S(7(6|3) 



tiuu + 11 + 4 



" [^uu + 2 + ^ujr] 91 + (1 + 92) (?^ r + 2 H 

-2i (T ss -2T us + T uu )(l- qi ))c 2 



A 

sjr 



8i 2 
TV 3 



/C ujr (£) + 3F + 2Fgi) + 2/C sjr (£> - F) + 2K, UU F (2 + q x ) 
+IC US (D + F)(2 + qi ) + (1 + g 2 )(F - D) 



-2i ^ ss (£)-F) + 25 u „F-5 us (D + F) ) (2 + gi) 
^W„ jr f (D - F)(D + 3F) + 2 (D 2 + F 2 ) Ql 



+ jH S3r (D-F) 2 + l -H uu [ IF' 2 +(D - F)(D+ F) qi 



--7i us ( 4p 2 -DF + F 2 ) + (D 2 + ADF - F 2 ) q x 



+-(1 + q 2 )(D 2 + 2DF + 5F 2 )H ujr + 1(1 + g 2 )(£) - F) 2 7i: s 
i? ss (D - F) 2 - 4R US F(D-F)+ 4R UU F 2 ) ( Ql - 4) 
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+ 21 sjr (l + qi ) + (1 + g 2 ) + 2I S 



+^u jr i(D-F)(D + 3F) - 4F 2 q x 



+ l -H sjr (D-Ff-^n uu {D' r F*-2F^, 
-jH ua (d 2 - 2DF -F 2 + 2(D - F)F Ql \ + |(1 + q 2 )(D - F) 2 (h ujt + 2H sjr ) 



1 



-- Rss (D - F) - AR US F{D-F)+ AR UU F 1 (1 + 2gi) 



(C33) 



U s -v {lMt . . . v^S^yUs- ((5 + A qi )Aa ( n a) - 2(-l + qi )A& 



+ - 



5„. _l£fe- 



6 

5i 

81 



X (l+W§(6|3)) 



(a) 



WA + W? a (2 + gx) + W£ (13 + 11 «i) 



+ 



+4 W£ r ] (1 + gi ) - (1 + q 2 ) (2 H% r + H% r ) - 2i (T ss -2T US + T uu ) (2 + q^C 2 
2K ujr {D -F)(l + qi) + K sjr {D + F + (D + 3F) qi ) + 2K ss F(l + 2 9l ) 



Si 2 



9 V 3 



+JC US (D + F)(l + 2 Ql ) + (1 + q 2 )(D - F) (2K ujr + t s 



-2% \S UU (D - F) + 2 S ss F - S us (D + F)j (1 + 2 gi) CAc { n a) 
+Aal a) ^ - J ujr gi + 5 J sjr (l + gi) + (1 + q 2 ) 

~H ujr (D - F) 2 (1 + qi ) - % -H sjr (d 2 - 2DF + 5F 2 - (D - F)(D + 3F) Ql 



+-7i as i^D 2 - 5 F 2 - 4 F 2 qi j + -H us {(3D - 5F)(D -F)+4(D 2 -DF + F 2 ) Ql 
+ + q 2 )(D - F) 2 n uir + kl + q2 ){D 2 + 2DF + 5F 2 )H sjr 



- 4R SS F 2 - AR US F(D-F) + R uu (D - F) 2 (5 + 4 Ql ) 



- - 21 ujr q 1 + l s j r (\ + gi) + (1 + q 2 ) 



-3 n ujr (D — F) (1 +qi) — - H sjr [{D + F) + {D-F)(D + 3F) qi 

+ | H ss (d 2 -F 2 + (D 2 + F 2 ) q^j + 1(1 + q 2 )(D - F) 2 (2 W UJ > + H sjr ) 
2i 



H us D 2 - 4DF + F + (D 2 - 2DF - F ) gi 



-- ( 4i? ss F - AR US F(D-F)+ R uu (D - F) 2 ) (1 - gi ) 



(C34) 



+ 



4V3 

Ux°v{y. ■ ■ ■ ^-i^„}^a°(-1 + gi) 



(a«w - 2A/3W) x [1 + \ (w£j m) + w. 



SU(6\3)J 



V3.f 2 



54 7 " 5 



Oi_/A 1 n_/A 1 n_/A \ s~i2 



27 



9 V 3 



2IC US (2D + 3F) - K, uu (D - 3F) + 2K ujr D + K sjr (D + 3F) 



i[2S uu (2D-3F) + S US (-5D + 3F) + S SS (D + 3F) 



CAc^ 



+Aai o) -l UJ r - ^H ujr (D - 3F)(3D - F) + ^H sjr (D - F)(D + 3F) 
+ fiuuD(D + F)+ l -H us D(D -F) + ^ (r ss {D - F) (D + 3F) 
-4R US (D 2 - 2DF + 3F 2 ) - AR UU F {-2D + 3F)) 



- X -T uir - l -H ujr (D - 3F)(D + F)- % -H sjr (D - F)(D + 3F) 



+ ^H UU D(D - F) + jH us DF - ± (r ss (D - F) (D + 3F) 
-AR US (D 2 - 2DF + 3F 2 ) - AR UU F {-2D + 3F) 



)]}■ 



(C35) 



4. Isoscalar unpolarised matrix elements 



(p\0^ n \p) = U p v, Q . . . «„„ tf p (a« + /#>) x (l + (1 - 5 n0 )W^ {m ) 



• ro 1 * I In 



+ 



P 



z n uu "+" n ujr I u 



(C36) 



i H ujr (5 D 2 - 6 D F + 9 F 2 ) + 4 iU uu D (D - 3F) - 3R UU (D - 3F) 2 



Y 



(A \O ( °l.., n |A°) = U A ov^ . . . v^U A o (a« + /?«) x (l + (1 - 6 n0 ) W$ (6|3) ) 



UaoV^o ■ ■ .V^UAojl - S na ) j . | (s) 



-ryA i a/A , ay A \ ,^2 



+ 



jH ujr (7D 2 - \2DF + 9F 2 ) - l -H sjr (D + 3F) 2 



P 

2? 27 
+ y W«„ (# 2 - Y2DF + 9F 2 ) + -W„ S (5C 2 - 6DF - 9F 2 ) 

-URss (D + 3F) 2 - 4R US (2D - 3F)(D + 3F) + 4R UU (2D - 3F) : 



(C37) 



(S + |0^... MB |S+) = Uv+vn . . . v^U^ (a« + /#>) x (l + (1 - 5 n0 ) W^ (6|3) ) 



+ 



Ux+v^ . . .^„L/ S +(1 - (5 ra0 ) / i 
P 



Tn 



is) 



^"ilu ^ ^~^us ^ ^sjr ^ujr 



-2 1 (T SS -2T US + T UU ) )C 



+ 



2 i H ujr (D 2 + 3F 2 )-3i H sjr (D-F) 2 + 2i H UU (D 2 - 3F 2 ) 



+2 i U US (D 2 - 6DF + 3F 2 ) - 3 ( R ss (D - F) 2 - 4R US F (D - F) + AR UU F 



)]} 
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(C38) 



(E-\0^ n \E~) = U B -v ll0 ...v lln U s -^+^) x (l + (1 - S n0 ) W^ (6|3) ) 



+ 



Us-Vfia ■ ■ -V^Us-il - 5 n0 ) I i I , . a 



+ 



-2i {T ss -2T us +T uu )^jC 2 

3 i Hu 3 r (D - Ff - 2 i Usjr (D 2 + 3F 2 ) + 2 i W SS (L> 2 - 3F 2 ) 

+2 i H us (i^ 2 - 6DF + 3F 2 ) - 3 (aR ss F 2 - 4R US F (D - F) + R uu (D - F) 2 ^j j . 

(C39) 



(«i s) +/3i s) ) 



{E^OW.^JA ) = 0. 



(C40) 



5. Isoscalar helicity matrix elements 



U p v { t l0 ...v, n _ 1 S lin} U p (Aatf + A(3M) x (l + W$ (6|3) ) 



/ 2 1 9 \ /n 5 



+ 



(Ac^+A/^)) 



L n uu + rt Mjr I ° 



4? 



W ujr (5 D — 6 D F + 9 F ) — -H uu D(D — 3F) + R UU {D - 3F)< 



(A \d%..„ n \A ) 



U A oV {flo . ■■Vp n _ 1 S lin yU A o (a4 s) + a/?w)x(i + >v^ ) (6| 3 ) ) 

/ A7 ( S) _ A^\ ^ +<> ^ 2 



(AaW + A/3( ; 



jH ujr {ID 2 - 12DF + 9F 2 ) + l -H sjr (D + 3F) 2 
~H UU (D 2 - 12DF + 9F 2 ) - ^H US (5D 2 - 6DF - 9F 2 ) 



+ - ( R ss (D + 3F) 2 - AR US {2D - 3F){D + 3F) + AR UU {2D - 3F)< 



(£+|0<E>...,JE+) = U^v^.-.v^S^yU^^Aa^+A^) x (l + wg? (6|3) ) 



[27 



») 



Act, 



+ 



x + 5 7Y^ S + 2 7Y^- r + - 2 i (T ss - 2 T„ s + T utt ) ^ C 



y W ujr (£ 2 + 3F 2 ) + iU sjr {D -Ff-j H UU {D 2 - 3F 2 ) 



(C41) 



(C42) 
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2i 



H US (D 2 - 6DF + 3F 2 ) + R ss (D - F) 2 - 4R US F (D - F) + 4R UU F 2 



)]}' 

(C43) 



+ (Aa«+A/#>) 



^2 + + + 5 H^s — 2i (T" ss - 2 T u;s + T utt ) ^ C 2 
i W„ jr (D - F) 2 + | W sjr p 2 + 3^ 2 ) - | H SS (D 2 - 3F 2 ) 



- % H us (D 2 - 6DF + 3F 2 ) + ( 4R SS F 2 - 4R US F (D - F) + R uu (D - Ff 



p \d^ n \A ) = 0. 



(C44) 
(C45) 



6. Trans versity matrix elements 



(p\o T „ na \ P ) 



U p v {» ■■ - v [^}S a ]U p ((5 + yi)5a n + 2(1 + 2y 1 )5(3 n) ^ ( l + w ip) 



G 



(i + W&s)) 



+ 



U p v^ . . .v^ n yS a ]U p f hi 



+ 



+5a, 



f 2 

8i [2 
9"V3 
1 



81 ( 5ln 



3 Sa r 



SJ7(6|3), 

3Ht+n* r ) (5 + yi )C 2 



KZ UU (D + 3F)(-1 + yi) - /C„ jr (3£> — F — 2F Vl ) 



C8c n 



G 



(5 + yi) J, 



-ujr 



2i( 



+ -H ujr [3D z ~6DF + 5F z + (D 2 + F 2 ) yi 



M 1 
--H uu D(D-2F- Fy x ) + -R UU (D - 3F) 2 (5 + yi ) 



+S(3 n 



-(1 + 2yi) ( T<y> - 2 / (';„ 



jH ujr (D 2 + F 2 (\ +2. Vjl ) 



~H UU D(D-2F- Fy t ) + ^R UU (D - 3F) 2 (1 + 2^ 



(C46) 



(A°I^,..^|A°) - £/A °^---^f" l[/AO(1 + i/l) ((l + 2^ + 2^) x (l + W, 

, UAQV^g ■ ■ ■V[ P , n }S a ]UAo( 1 + Vi) f 5i / 3 Joy 

+ / 2 \54 5 



(A) \ 
SC/(6|3) I 



(1 + y) + Ut (1 + 4y) + H% r (1 + 2y) j C 2 



■-a/- 

9 V 3 
1 



- IC US (D - 3F) + 2IC UU D) (1 — 2y) — - 3.F + 4Dy) 



+ 21 sjr y - 2iQ uu - 4iQ ss y 
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+z^H uu [ D 2 + \2DF - 9F 2 - 6(D - 3F)(D - F)y 



~H US [ (D - 3F)(5D + 3F) + 2(5D 2 - 9F 2 )y 



+ — (R ss (D + 3F) 2 - AR US (2D - 3F)(D + 3F) + AR UU (2D - 3F) 2 j (1 + 2y) 
+6f3 n - X - (l ujr - 2 1 Q u ^j + ^H ujr ({D - 3F) 2 + 8D 2 y^j + ^H sjr (D + 3F) 2 
-\Uuu U(D - 3F)(D — F) — AD 2 y) + l -H us (d 2 + 9F 2 - 12D(D - F)y 



+ — [R ss (D + 3F) 2 - AR US (2D - 3F)(D + 3F) + AR UU (2D - 3F) 2 J j j . (C47) 

Ut,+ V{^ ■ ■ ■ v[» n }S a ]Uv+ ((5 + y 2 )Sa n + 2(1 + 2y 2 )6(3 n ) / (s) \ 

g x y 1 + yv su(e\3)J 

+ p 1 si 1^7" — J I H uu ( 2 + Wz) + W us ( 13 + 2 2/2) 

+4 U% r + U% r (1 + y 2 ) -2i(T ss -2 T us + T uu ) (2 + y 2 ) S j C 2 

~?Vl Ku ^ D + F ~ 2Fy2 ^ + 2IC ^(D — F) + UlC uu F + JC US (D + F)\ (1 - y 2 ) 



+5a Tl 



-2i(S ss (D-F) + 2S UU F-S US (D + F) j (1 - y 2 ) 

g lujr — 10i Q MU + \lsjr — 2i Q S s\ y^j 

~H ujr (d 2 - 2DF + 5F 2 + 2(D 2 + F 2 )y 2 ) + |w sjr (D - F) 2 



C<5c„ 



-gW™ ( F 2 - 5F 2 + (D-F)(D + F)y 2 



-^H us ( 3£ 2 - 8DF + 5F 2 - (D 2 + ADF - F 2 )y 2 



+ - ( R ss (D - F) - AR US F(D-F)+ AR UU F 2 ) (5 + y 2 ) 



+S(3 n 



.'72 



(D + F) 2 + AF 2 y 2 



2i 



--U uu [D 2 -F 2 -2F 2 y 2 



2i 



H us D 2 - ADF + F — 2F(D - F)y 2 



i„. „, 2 , 1 



+-n sjr (d - f) 



R ss (D - F) 2 - AR US F(D-F)+ AR UU F 2 )(l + 2y 2 ) 



(C48) 



Ue-V{^o ■ ■ ■ v [^} s a ]U s - ((Vi + hy 2 )Sa n + (A Vl + 2y 2 )5f3 n ) 



(i + wf^eis)) 



+ 



f 2 



6 

| Uln-^) (<>(2/i+y2)+4< >y2 +Wf s ( /yj ; 2,/o) 



+W^ S (2 yi + 13y 2 ) - 2z (T ss -2T US + T uu ) ( Vl + 2y 2 ) 



C 2 



8i 2 

TV 3 



2/C ujr ( J D - F)y 2 - /C sjr 2F Vl - (D + F)y 2 - K US (D + F)( Vl - y 2 ) 



31 



~2JC ss F{ yi - y 2 ) + 2i S uu {D-F) + 2S SS F- S us (D + F) ( m - y 2 ) 



CSc n 



-5a n 



- g ( [lujr - 2 i Quu] 2/1+5 [l s j r -'lit 



in 



~H ujr (D - F f y 2 + l -H sjr ( 2(D 2 + F 2 ) yi + (D 2 - 2DF + 5F 2 )u> 



-H ss ( (D - F)(D + F) Vl + (D 2 - 5F 2 )y 2 



--Hus ( (D 2 + ADF - F 2 )yi - (3D - 5F)(D - F)y 2 



( , , AR SS F - AR US F{D-F)+ R uu (D - F) ) (tfi + 5y 2 ) 
i ( 2 [l ujr - 2 i Q uu ] yi + [l s j r - 2 i Q ss ] y 2 J + ^ H U j r [D - F) 2 y 2 



3 • ?r 



AF 2 Vl + (D + F) 2 y 2 



2% 



H ss [2F 2 y 1 - (D 2 - F 2 )y 2 



2i 

T 
i 



Tins ( 2F{D - F) yi - (D 2 - ADF + F 2 )y : 



AR SS F - AR US F(D-F) + R uu (D - F) (2 Vl + y 2 ) 



(C49) 



(S°|^ ...^JA° 



+ 



4^3 

UsoV{no ■ ■ ■V[ t 2„}S a ]UAo(-l +m) ( _ 5i 
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(6a n - 25f3 n ) x 



1 +^H A a ) ( 6 |3)+W^ ) (6|3)) 



V3f 2 
Ai [2 

TV 3 



6ln 



3 8a n 



2K US (2D + 3F) - K uu (D - 3F) + 2K uyr D + tC s j r (D + 3F) 



i 2S UU {2D-3F)+S US (-5D + 3F) + S SS (D + 3F) 



+Sa n 



C5c n 



2 1 Quu)- -^Hujr (D - 3F)(3D - F) + —H sjr (D - F)(D + 3F) 



12 



12 



l -UuuD{D + F)+ l -Hu S D(D - F) + i( (I) - - F ) ( D + 3F) 



-4i? us (L> 2 - 2.DF + 3F 2 ) - 4i? UM F (-2F> + 3F) 



+6(3„ 



- gW UJ> (-D - 3F)(D + F) - -H sjr (D - F)(D + 3F) 



~H UU D(D — F) + %U US DF — i ( R ss (D — F) (D + 3F) 
3 3 6 



-AR US {D 2 - 2DF + 3F 2 ) - AR UU F {-2D + 3F) 



(C50) 



APPENDIX D: RESULTS IN SU(2|2) QxPT 

In SU(2|2) quenched xPT, the Lagrangian Eq. \12\ receives additional contributions since the theory has no axial 
anomaly and the singlet meson field remains light. Thus, 



{Bv -vb)-i [r v -vrJ+A [T% 



.32 



+2a (BS^BAfj) + 2/3 (BS" + 27i (t'VAX 

+2 7 (gS"B) str (A M ) + 2 7 ' (T"S"T tf ) str (A„) , (Dl) 

with two additional couplings 7 and 7'. There is no relation between the other couplings in Eq. (|D1|> and those in 
(PQ)xPT (though we use the same notation for convenience). 
Defining 



f a = diag(l, -l,q,-q) 



(D2) 



f =diag(l,l,l,l) 



(D3) 



and 

f T = diag(l,yi, 1/3,2/3) , (D4) 

the quenched twist-two operators correspond (for the most part) to those given in Eqs. I|18(l - (|20|l with the 

replacement r — * f everywhere. Again, the LECs (a4 etc) occurring in the quenched theory are different from 
those in SU(4|2) but with this caution, we use the same notation. For the transversity operators one needs additional 

operators proportional to str (^f^i J : 

[ Sa 'n v {no ■■■ v [fi„} {SS a] B) +S-/' n v {fl0 .-.%„} (T v S a] T v ) +8a' n v {lM) ■ ■■v Un _ 2 {T a n ^S [a f an]} )] str . (D5) 

With these definitions it is then easy to calculate the quenched matrix elements in the isospin limit. The wave- 
function renormalisation is 



Wsu(2\2) = -p ( (91 + 9 a) (67 - 31 + 2g A ) H uu + 2 g AN H uu + (gi + g A ) m H n 
The isovector, unpolarised matrix element is 



t \N) = ^U N v^...v Un U N ^2a ( ^ -(3 ( n a ^ x (l + (1 - 6 n0 )Wsu(2\2^\ 



(D6) 



(D7) 



+ JJ2 UnV vo ■ ■■^ n U N (l - 5 n0 )\4:ig AN ( 7 



(a) 



i(2aW-/9W) x 



QgAlT-Luu + (31 + 3a) w W,, 



+a 



(a) 



3i 







~3i 1 




+ A a) 


_-2 9l { 



and the isovector, helicity matrix element is 



(NpW..^ \N) = ^U N v {lM> . . . v^S^Un - A&A + W su(m ^j 



8i /2 . A _ . 20 .0 I a (n\ Act 
y Y 3 (3i + 4.g j4 )/C ull 3AArAc„ + y 1 3 AAr I A 7 ^ ' 



(a)' 



(D8) 



+- ( 2Aa£° - Af3^ ) x 



31 - 2(4 7 + 34) 
The isoscalar, unpolarised matrix element is 



63^7^™ + (31 + 9A) 2 mlH r f^ 

g 1 + 2{2 1 -g A ) H„ 



H uu - 2^31 



<JV|0^...^ |iV) = U N v U0 . . . v Un U N ( + ffi) ) x ( 1 + (1 - S n0 ) W sum) 



M 1 fl( s ) 



(D9) 
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+ 



OnOj 

P 



*(*) 



2 i 9an^-uu 



- i (a { n s) + /#°) (91 + 9a) ((67 - 51 + 2<m)W„„ + (fli + &0m§«„/, u1 ^ J • 



The isoscalar, helicity matrix clement is 



(D10) 



+ 



i7iyf{ w -..^„- 1 5 fa }^ f ( A7 ( S ) _ Act}* 



,/' 2 



1 



10 i 



■ 9an7~L II 



+ l - \AaM + A/J<°> ) ( 9l + g A ) ( (67 - 91 + 2g A )H uu + (.gi + <M)mX', 



Finally, the transversity matrix elements are 



UNV{y a ---V[^}S a] U. 



N 



+ 



UnV{^ ■ . .V[fj. n }S a ]U N 
P 



(5 + yi)Sa n + (2 + 4yi)S/3 n + 6(1 +yi-y 2 - m)5a' n ) ( 1 + W SU (2\2) 
(6 + 2 yi - y 3 - y 2 ) (5 Sj n - 3 5<r n ) 
+6 (1 + yi-ys-fc) (5*7;-3JO 



— — T-/ 

2y ^uu 



2 

5 a at 



4i /2 



9 V 3 

+5a n 



/C, 



4g A (l - Vl ) + 5l (4 + 2yi - 3y 2 - 3y 3 ) 



9 AN Sc n 



(-1 - 2/1 + 2/2 + 2/3)2"™ - -TOq(5 + 2/1) 2 "))' 



+ 12 ( -Si (2 + y 3 + 2/2) + 4 5A [,g A (2 + 2 & - y 3 - y 2 ) + (5 + 2/1) 7] 



+2. 9l [,g A (3 + Vl - y 3 - y 2 ) + 2 (5 + 2/1) 7]^ W„ 



+ 7^(.9i +5a) 2 ™o( 5 + 2/i) w >)'^ 



(-1 - 2/1 + 2/2 + yz)luu - g m o(! + 22/i)^',t, 



+ — ( 4 5l (2 7 + ,g A + 473/1) + 4.g A [2 7 + ffA + (47 + 9A)yi] 



+5a' 



+.g 2 [2yi - 3(2/2 + 2/3)] J W«u + -(.gi + .gA) 2 ^(l + 2yi) 
2 

4(2/2 + 2/3)2"«« - g^ol 1 + 2/i - 2/2 - h)l n >,uu 



■gCfli + Su) ( .gi - 2(37 + g^) ) (1 + 2/1 - 2/2 - 2/3) K 



+ g(ffi + g^) 2 ™ 2 ,^ + 2/1 - 2/2 - 2/3) Wr/ 



(Dll) 
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